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ABSTRACT 

Aims. We determine the solar abundance of phosphorus using CO^^BOLD 3D hydrodynamic model atmospheres. 
, Methods. High resolution, high signal-to-noise solar spectra of the Pi lines of Multiple! 1 at 1051-1068 nm are compared to line formation 
computations performed on a CO^BOLD solar model atmosphere. 

Results. We find A(P)=5.46 ± 0.04, in good agreement with previous analysis based on ID model atmospheres, due to the fact that the P i lines 
of Mult. 1 are little affected by 3D effects. We cannot confirm an earlier claim by other authors of a downward revision of the solar P abundance 
by 0. 1 dex employing a 3D model atmosphere. 
J Concerning other stars, we found modest (< 0. 1 dex) 3D abundance corrections for P among four F-dwarf model atmospheres of different 
metallicity, being largest at lowest metallicity. 

Conclusions. We conclude that 3D abundance corrections are generally rather small for the Pi lines studied in this work. They are marginally 
relevant for metal-poor stars, but may be neglected in the Sun. 
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1. Introduction 

Phosphorus is an element whose nucleosynthesis is not clearly 
understood and whose abundance is very poorly known out- 
side the solar system. It has a single stable isotope ^'P, and 
its most likely site of production are carbon and neon burning 
shells in the late stages of the evolution of massive stars, which 
end up as Type II SNe. The production mechanism is proba- 
bly through neutron ca pture, as it is for the parent nuclei ^"^Si 
and ^"Si. According to IWoosIey & Weaver I (Il995h there is no 



significant P production during the explosive phases. 

Like for other odd-Z elements (Na, Al) one expects that 
the P abundance should be proportional to the neutron excess 

As metallicity decreases the neutron excess decreases and 
these odd-Z elements should decrease more rapidly than neigh- 
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' as defined in lAmet 3 h97lll : T] = {n„ - np)/(nn + iip), where 
represents the total number of neutrons per gram, both free and bound 
in nuclei, and iiy, the corresponding number for protons 



bouring even-Z elements. One can therefore expect decreasing 
[""P/^^Si] or [""P/'^Mg] with decreasing metalHcity. Such a be- 
haviour is in fact observed for Na (Andrievsky et al. 2 0071) . 
however at low metallicity [Na/Mg] seems to reach a plateau 
at about [Na/Mg] ~ -0.5. One could expect a similar behaviour 
forP 

Phosphorus abundances in stars have so far been deter- 
mined for object classes not apt t o study the chemice tl evolu- 
tion: chemica llv pecuhar stars (e.g. Castelli et al.ll 19971: Fremat 
& Houziaux Il997 . for a review of observations in Hg-Mn 
stars see Takada-Hidai 1991 an d references therein), h orizontal 
branch star s tBehr et al. < J 9991: Bonifacio et al.lll995 LBaschek 



& Sargent 1 1976h. sub-dwarf B tvpe stars jOhl et al 



1982h. Wolf-Rave t stars dMarcohno et aljl2007h 



20001: 



Baschek et al.l , . 

P G 1159 stars d.Tahn et alJl2007h. and white dwarfs ( Chaver et 



al. l2005l : lDobbie et al. II2005I: IVennes et al. Ill996l) . Such stars 
are not useful to trace the chemical evolution because either 
they are chemically peculiar or they are evolved and may have 
changed their original chemical composition. A few measure- 



ments exist in O super-giants (Bouret et al. 20 U5t Crow ther et 
al. .2002) and B type stars (Tobin & KaufmanlTl 1983). Kato 
et al. (Il996h measured P in Procyon. ISion et al. I (Il997h in the 
dwarf nova VW Hydri. 

The abundance of P in F, G and K stars, which could be 
derived from the high excitation IR Pi lines of Mult. 1 at 
105 1-1068 nm has never been explored, due to the few high 
resolution near IR spectrographs available. With the advent of 
the CRIRES spectrograph at the VLT, the situation is likely to 
change and it will be possible to study the chemical evolution 
of P in the Galaxy using these long-lived stars as tracers. In this 
perspective it is interesting, as a reference, to revisit the solar P 
abundance in the light of the recent advances of 3D hydrody- 
namical model atmospheres. There are several determinations 
of the solar P abundance in the literature. Given the difficulty 
of this observation, it is not surp rising tha t P was ab sent from 
the seminal work of lRussell I (11929 ). Lambert & Warner (1968) 
measured an abundance A(Pj3 - 5.43, using 13 hnes of Pi, 
of which four belong to Mult. 1 . The o scillator strength they 



lable 1. Infra-red phosphorus lines considered in this work. 



used f or these lines are ver y close to the | Biemont & Grevesse . 
I973h values (see Table [3Tl. lLambert & Luckl (Il978h measured 



A(P) = 5.45+0.03, which is essent ially the photospheri c abun- 
dance adopted in the authoritative lAnders & Grevesse (Il989l) 
compilation (A(P)= 5.45 + 0.04). Based on their ab initio com- 
puted oscillator strengths Biemont et al. ( I994I) derived A(P) 
= 5.45 + 0.06. With computed oscillato r strengths, corrected to 
match measured energy level lifetimes, iBerzinsh et al. I (Il997h 
deriv ed 5.49 + 0.04. In their compilation Grevesse & Sauval 
IT998) adopt A(P) = 5.45 + 0.04. As can be seen, all de- 
terminations of the solar photospheric phosphorus abundance 
are in agreement, withi n the stated u ncertainties. In disagree- 
ment with earlier work, lAsplund et al. (2005) obtained a sig- 
nificantly lower value of A(P)= 5.36 ± 0.04 using a 3D solar 
model atmosphere computed with the Stein & Nordlund code 
(see Nordlund & Stein 1997). This result might be immediately 
interpreted as indication that 3D abundance corrections lead to 
a lower solar phosphorus abundance by ~ 0.1 dex, but other 
factors, such as the choice of logg/ values, can also explain 
this low A(P). The A(P) determinations found in literature are 
summarised in Table |2] 

In the present paper we use a solar 3D CO^BOLD model to 
re-assess the solar phosphorus abundance. For comparison we 
also consider ID solar models. Beyond the Sun, we further dis- 
cuss abundance corrections obtained for a number of F-dwarf 
models of varying metallicity. 



2. Atomic data 

We consider five IR Pi lines of Mult. 1 (see Table [U. Several 
determinations of the oscillator strengths for this multiplet are 
available in the literature (see Table [3]). Generally, all values 
are rather close, differences are smaller than 0. 1 dex for any 
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line except the 1051.2nm line for which the maximum differ- 
en ce is 0.18 dex. In this w ork we adopted the values computed 
by IBerzinsh et al. I (Il997h . corrected semi-empirically, which, 
in our opinion, are the best currently available. Had we cho- 
sen another set of val ues, with the notable exception of the 



Berzinsh et al. (Il997 l) ab initio values, the difference in the 



derived phosphorus abundance would be 0.05 dex at most. 



3. Models 

Our analysis is based on a 3D model atmosp here computed 
wi th the CO^BOLD code dFrevtag et al. II2OO2I: Wedemever et 
al. |2004 . In addition to the CO^BOLD hydrodynamical sim- 
ulation we used several ID models. More details about the em- 



ployed models can be found inlCaffau et al.l (l2007i) and Caffau 
& Ludwig (2007). In the following we will refer to the ID 
model derived by a horizontal and temporal averaging of the 
3D CO^BOLD model as the (3D) model. The reference ID 
model for the computation of the 3D abundance corrections 
is a related (same T^s, logg, and chemical composition) hy- 
drostatic ID model atmosphere computed with the LHD code, 
hereafter IDlhd model. LHD is a Lagrangian ID hydrodynam- 
ical model atmosphere code, which employs the same micro- 
physics as CO^'BOLD. Convection is treated within the mixing- 
length approach. In the LHD model computation, as in any ID 
theoretical atmosphere model, at least one free parameter de- 
scribing the treatment of convection must be decided upon: we 
use a mixing-length parameter of a=1.25 in the solar case, a 
value of a= 1 . 8 in [M/H] =-2.0, and g^I.Oin [M/H] =-3 .0 mod- 
els. We employ the formulation o f lMihalasl(ll978h for mixing- 
length theory, and we neglect turbulent pressure in the momen- 
tum equation. We also considered the solar ATLAS 9 model 
computed by Fiorell a Castelli H and the empirical Holweger- 



Miille r solar model (iHolweger 
1974. 



I967t IHolweger & Miiller 



The spectrum synthesi s codes employed are Linfor3E0for 
all mod els, and SYNTHE (lKuruczlll993bl.l2005al) . in its Linux 
version dSbordone et al. Il2004l:lsbordone ll2'o05h . for the fitting 
procedure with ID or {3D) models. The advantage of using 



A(P) = log(A'p/A'H) + 12 



' http://wwwuser.oats.inaf.it/castelli/sun/ap00t5777g44377klasp.dat 
http://www.aip.de/ mst/Linfor3D/linfor_3D_manual.pdf 



SYN IHb, with respect to the present version or LinrorJU, is 
that it can handle easily a large number of lines, thus allowing 
to take into account numerous weak blends. 

The 3D CO^BOLD solar model we used is the same we de- 



scribe in lCaffau & Ludwig I (l2007h . It covers a time interval of 
6000 s, represented by 25 snapshots; this covers about 10 con- 
vective turn-over time scales, and 20 periods of the 5 minute os- 
cillations which are present in the 3D model as acoustic modes 

of the computational box. 

Like in our study of sulphur (seel Caffau & Ludwig 112007 ). 
we define as "3D abundance correction" the difference in the 
abundance derived from the full 3D model and the related 
IDlhd model, in the sense A(3D) - A(IDlhd), both synthesised 
with Linfor3D. We consider also the difference in the abun- 
dance derived from the 3D model and from the (3D) model. 
Since the 3D and (3D) models have, by construction, the same 
mean temperature structure, this allows to single out the effects 
due to the horizontal temperature fluctuations. 

4. Data 

As observational data we used two high-resolution, high signal- 
to-noise ratio spectra o f disk-centre solar intensity: that of 
Neckel & Labs I d 19841) (hereafter refeiTed to as the "Neckel 
intensity spectrum") and that of Delbouille, Roland, Brault, & 
Testerman (1981) (hereafter referred to as the "Delbouille in- 
tensity spectrum" |£|)- W e also used the solar flux spec trum of 



Neckel & Labs I (119841) and the solar flux spectrum of lKurucz 
( 2005H4 referred to as the "Kurucz flux spectrum". 



5. Data analysis 

We measured the equivalent width (EW) of the lines using the 
IRAF task splot. We then computed the phosphorus abun- 
dance for the solar photosphere from the curve of growth of 
the line in question calculated with Linfor3D. As a confirma- 
tion of our results we fitted all observed line profiles with syn- 
thetic profiles. We found a good agreement in comparison to 
the abundances derived from the EWs. The remaining differ- 
ences we found are probably related to the difference in the 
continuum opacities used by S YNTHE with respect to the ones 
used by Linfor3D. 

The line profile fitti ng was done using a code, described 



in iCaffau et al. I (|2005l) . which performs a x minimisation 
of the deviation between synthetic profiles and the observed 
spectrum. Figure [1] shows the fit of the solar intensity spec- 
trum obtained from a grid of synthetic spectra synthesised with 
Linfor3D and based on the 3D CO^BOLD model. The close 
agreement between observed and synthetic line profile indi- 
cates that the non-thermal (turbulent) line broadening is very 
well represented by the hydrodynamical velocity field of the 
CO'^BOLD simulation. 
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Fig. 1. The observed disk-centre Neckel intensity spectrum 
(solid line) is plotted over the fit (crosses) obtained using a grid 
of synthetic spectra based on the CO^BOLD model computed 
with Linfor3D. 



6. Results and discussion 

6.1. The solar P abundance 

The different sets of log gf values for the selected P i lines 
are very close, so that the derived solar A(P) is quite insen- 
sitive to this choice. The solar phosphorus abundance varies 
by 0. 05 dex, from th e high est value obtained with the data 
set of Biemont et al.l ( 1994) to the lower one obtained using 
Biemont & Grevesse lll973i). the corrected value of Berzins h 
et al. (Il997ir or the one from lKurucz & PevtremannI ([19751) • 
A lower value of the solar pho sphorus abundance (0. 09 dex be- 
low the one obtained using the lBiemont et al.l (119941) oscillator 
strength) can b e obta ined using the non-corrected values from 



Berzinsh et aU ( [l997h . 



With the (corrected) \og gf of iBerzinsh et al. I ( 119971) we 



obtain A(P)=5.443 + 0.058 for the flux spectra considering the 
whole sample of five lines, A(P)=5 .426+0.064 for the intensity 
spectra at disk-centre. The line at 1068. Inm yields an abun- 
dance which is 1 .6 cr below the mean value, while the other four 
lines are within one cr. If we remove this line from the compu- 
tation, the standard deviation drops by a factor of two and we 
obtain A(P)=5.467 + 0.029, 5.450 + 0.042 for flux and inten- 
sity spectra, respectively. The differences in the A(P) determi- 
nation from intensity and flux spectra could be due to NLTE 
corrections which usually are different in intensity and flux. 
Unfortunately, there are no NLTE computations available for 
these lines. Our recommended solar photospheric abundance is 
A(P)=5.458 + 0.036, very close to all previous ID determina- 
tions, however 0. 1 dex higher than the value found bv Asplund 
et al. ( 2005) using a 3D model. 



^ http://bass2000.obspm.fr/solaru_spect.php 
* http://kurucz.harvard.edu/sun.html 



The result of lAsplund et al.l (l2005h is difficult to interpret 
since, according to our m odel, the 3D co r rection s for the Pi 
lines in the Sun are small. lAsplund et al.l (120051) do not pro- 
vide detailed information on the lines used or of t heir a dopted 
EWs. They adopted the \og gf of IBerzinsh et al. I ( Il997h . how- 
ever they do not specify if they used the "corrected" or ab initio 
values. If they used the ab initio values the difference is easily 
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Fig. 2. Relative horizontal RMS temperature fluctuations (on 
surfaces of equal Rosseland optical depth) in the solar 
CO^BOLD model. 



understood, otherwise we must conclude that the difference is 
due either to a difference between their 3D simulation and the 
CO^BOLD one, or due to a difference in the adopted EWs. We 
point out here that our measured EWs are in quite good agree- 
ment withbothjhose of La mbert & Warner ( 196 8). Lambert & 
Luck d 19781) and with fliose of lBiemont et all 11994 

A detailed summary of the solar phosphorus abundances 
we derived from individual lines, using different model atmo- 
spheres, in given in Table HI Since the considered lines are not 
truly weak, the 3D corrections are unfortunately sensitive to the 
adopted micro-turbulence, and hence must be interpreted with 
care. The 3D corrections due to horizontal temperature fluctu- 
ations only, indicated by the 3D-(3D) difference, are slightly 
positive throughout. In all lines, the 3D-(3D) difference for in- 
tensity spectra is systematically higher than for flux spectra, by 
roughly 0.02 dex. This is presumably due to the fact that the 
intensity spectra originate from somewhat deeper layers where 
horizontal fluctuations are larger. We note that the phosphorus 
lines are formed mostly in the range -0.5 < logjg Tioss < 0, 
where the horizontal temperature fluctuations increase with in- 
creasing optical depth (see Fig.|2]i. 

The 3D-(3D) difference increases with the excitation po- 
tential of the line's lower level. In fact, the 3D-(3D) correction 
is rather small, not exceeding +0.030 dex for f=LO km s \ and 
H-0.045dex for ^=1.5 kms ', which we consider as an upper 
limit. 

The 3D-1Dlhd corrections are systematically larger than 
the 3D-(3D) difference (by +0.01 to 0.02 dex), but of the same 
order of magnitude. From that we can conclude that the solar 
(3D) and IDlhd temperature profiles are very similar 



6.2. 3D effects on P for other stars. 

As shown above, we obtain only mild 3D effects in our solar 
phosphorus abundance analysis. They are comparable to the 
standard deviation (as listed in col. (9) of Table [J]) for the weak 



fines, and a factor tour larger tfian tfie standard deviation tor 
tfie strongest line, and fience tfie 3D effects are fiardly relevant. 
However, 3D effects sfiould become stronger for metal-poor 
stars wfiere tfie temperature profiles differ more strongly be- 
twe en (3D) and ID i .Hn models tfian for solar metallicity mod- 
els (lAsplundll2005l) . Very little work fias been dedicated to 3D 
effects on tfie pfiospfiorus abundance, and none for metal-poor 
stellar models. To cfieck tfie befiaviour of tfie 3D abundance 
corrections, we tfieoretically investigated for a few available 
3D CO^BOLD models tfie pfiospfiorus 3D- IDlhd relation in 
tfie resulting syntfietic flux spectra. 

We find, for tfiree metal-poor models (6300/4.5/-2.cO; 
5900/4.5/-3.0; 6500/4.5/-3.0), tfiat tfie difference in tfie mean 
temperature profile is tfie only factor wfiicfi contributes to tfie 
3D correction. In tfie metal-poor models, as in tfie solar model, 
tfie maximum contribution to tfie equivalent widtfis of tfie P i 
lines comes from a well defined region close to tfie continuum 
forming layers. Even tfiougfi tfie optical deptfi range contribut- 
ing significantly is somewfiat broader in tfie metal-poor models 
tfian in tfie solar model, tfie (3D) and IDlhd models are usu- 
ally ratfier close in temperature in tfiis region, and so tfie effects 
of tfie different mean temperature structure are not very pro- 
nounced for tfie considered Pi lines. Tfie 3D correction related 
to tfie difference in tfie mean temperature structure of tfie 3D 
and IDlhd model amounts to ^ +0.05 dex for tfie -2.0 metal- 
licity model, and to x +0.1 dex for tfie most metal-poor model. 

For a 3D CO^BOLD model of Procyon (6500/4.0/0.0) tfie 
difference in tfie temperature profile of tfie (3D) and IDlhd 
model contributes to more tfian fialf of tfie total 3D correction. 
In Procyon, Pi lines are formed at sfiallower optical deptfis 
tfian in tfie Sun by about Alogr w 0.3. In tfiis range tfie 
IDlhd model is cooler tfian tfie 3D CO^BOLD one. Tfie re- 
lated abundance correction depends on tfie excitation poten- 
tial (strengtfi) of tfie line, and for tfie strongest line it is still 
less tfian +0.030 dex. Tfie contribution to tfie 3D correction 
related to tfie fiorizontal temperature fluctuations ranges from 
+0.022 to +0.045 dex, increasing witfi tfie excitation potential 
(strengtfi) of tfie line. In tfie range wfiere tfie line is formed, tfie 
fiorizontal RMS temperature fluctuations in Procyon are of tfie 
order of 5% to be compared to 1% in tfie metal-poor models. 
Tfie total 3D correction is tfius less tfian +0.07 dex for aU fines. 



7. Conclusions 

Using tfie CO^BOLD solar granulation model, we fiave de- 
termined tfie solar 3D LTE pfiospfiorus abundance to be 
A(P)=5.46 +0.04, wfiicfi compares very well to all previous de- 
terminations obtained using ID models. Tfiis can be explained 
by tfie fact tfiat tfie P i lines of Mult. 1 appear to be ratfier in- 
sensitive to granulation effects, at least in tfie parameter inter- 
val explored by us (F and G dwarfs of metallicity from solar to 
-3.0). Tfiese lines can tfius be useful to investigate tfie cfiem- 



^ indicating r,ff/logg/[M/H] 



ical evolution or phosphorus in the Oalactic disk and in mod- 
erately metal-poor environments. Below [P/H]=-1.0, however, 
the lines have EWs smaller than 0.5 pm, becoming exceedingly 
difficult to observe. 

For the P i lines studied in this work, we find small (Sun) 
to moderate (Procyon, metal-poor stars) 3D abundance correc- 
tions. The sign of the corrections is found to be positive in all 
cases, meaning that the abundances resulting from a 3D analy- 
sis are larger than those obtained from a ID model. Note that 
this behaviour is consistent with the results obtained bv Steffen 
& Holweger C200Z) for S i, which has almost the same ionisa- 
tion potential as Pi (see their Fig. 8). In general, however, the 
magnitude and sign of the 3D effects depend on the properties 
of the absorbing ion (ionisation potential, the excitation poten- 
tial of the lower level, line strength and wavelength), and on the 
thermodynamic structure of the stellar atmosphere. The results 
obtained in this work for P i must therefore not be generalised 
to other elements and other types of stars. 
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lable 2. Spectroscopic determinations or the solar photo- 
spheric phosphorus abundance. 



A(P) 


O" 


Reference 


5.43 




Lambert & Warner (1968) 


5.45 


0.03 


Lambert & Luck ("1978) 


5.45 


0.04 


Anders & Grevesse ('1989) 


5.45 


0.06 


Biemont et al. (1994) 


5.49 


0.04 


Berzinsh et al. (1997) 


5.45 


0.04 


Grevesse & Sauval (1998) 


5.36 


0.04 


AsDlund et al. (2005) 


5.46 


0.04 


this work 



Table 3. Oscillator strength for the IR phosphorus lines (wave- 
lengths in nm) available in the literature. 



1051.2 


1053.0 


1058.2 


1059.7 


1068.1 


Reference 


-0.21 


0.20 


0.48 


-0.21 


-0.10 


Kurucz & Pevtremann (1975) 


-0.23 


0.18 




-0.24 


-0.12 


NIST 


-0.06 


0.31 


0.52 


-0.14 


-0.12 


Berzinsh et al. (1997) 


-0.13 


0.24 


0.45 


-0.21 


-0.19 


Berzinsh et al. (1997) corr. 


-0.10 


0.27 


0.49 


-0.17 


-0.15 


Biemont etal. (1994) 


-0.24 


0.17 


0.45 


-0.24 


-0.13 


Biemont & Grevesse (1973) 



(NIST: |http://physicsjiIsrgov/PhysRefData/ADS. I 



lable 4. Solar phosphorus abundances rrom the various observed spectra, adopting the (corrected) ioggj or IBerzinsh et al. I 
d 19971) . 



Spectrum 


Wave 


EW 




A(P) from 


EW 




a 


3D-1Dlhd 


3D-<3D> 




[nm] 


[pm] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 
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(11) 


( 1 7^ 
(IZ) 


(13) 


KF 


1051.2 


0.674 


5.482 




S Al\ 


J .H-OU 


S A7S 


0.015 


0.015 


0.019 


0.005 


0.010 




1 1 9 

iUJ 1 .z. 


u. / 




J.J I 1 


J.J iZ. 


J.JZ.L 


S S 1 A 

J. J lO 


U.Ul J 


(1 ill ^ 


(1 (17(1 

u.uzu 


U.UUj 


(1 (11 (1 
U.UIU 


NI 


1051.2 


0.811 


5.502 


5.476 


5.473 


5.498 


5.493 


0.012 


U.UZo 


(1 (1^7 
U.UjZ 




n n7A 

U.UZO 


DI 


1051.2 


0.822 


5.508 


5.483 


5.479 


5.504 


5.500 


0.012 


U.UZb 


(1 (n7 
U.UjZ 


n n77 
U.UZZ 


(1 (177 
U.UZ / 


KF 


1053.0 


1.330 


5.431 




41 9 


477 


S 41 7 


0.010 


0.019 


0.028 


0.005 


0.015 


NF 


LyJJJ.KJ 








S 41 ^1 


S 4^ 1 


S 49 1 


U.UIU 


n ni Q 


(1 (17 Q 

u.uzs 


U.UUj 


(1 (11 '^ 

U.Ul J 


NI 


1053.0 


1.520 


5.425 


5.392 


5.382 


5.417 


5.408 


0.010 


U.Ujju 


(1 (1/1/1 


n (17 
U.UZO 


(1 (1^/1 


DI 


1053.0 


1.550 


5.436 


5.402 


5.392 


5.427 


5.418 


0.010 




(1 C\AA 


n (17 s 

U.UZJ 


(1 (l^A 


KF 


1058.2 


2.280 


5.480 




AA9. 


'\ All 
J.^ 1 J 




0.012 


0.027 


0.042 


0.006 


0.023 


NF 


1 9 

iUJO.Z, 


9 1 9n 




S Al 1 

J .H-Z. 1 


S AC\f\ 

J .H-UU 


J .H-Zy 
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U.Ul J 


U.Uzj 


(1 (ITQ 


U.UUj 


(1 (17 1 
U.UZl 


NI 


1058.2 


2.330 


5.407 


5.363 


5.346 


5.393 


5.379 


0.012 
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(1 (1'^Q 
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n (170 
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(1 (1/1 T 


DI 


1058.2 
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1059.7 
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0.005 


0.010 
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(1 (11 (1 
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1059.7 
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5.431 


5.427 
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5.448 


0.012 
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(1 (ITd 

U.UjU 
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U.UZZ 


(1 (17A 
U.UZO 


DI 


1059.7 


0.751 


5.458 


5.433 


5.429 


5.454 


5.450 


0.012 


n (177 
U.UZ / 


(1 (IT 1 
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U.UZZ 
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1068.1 
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5.370 


5.364 


5.360 


5.368 


5.364 


0.014 


0.015 


0.019 


0.005 


0.010 


NF 


1068.1 


0.617 


5.330 


5.324 


5.320 


5.328 


5.324 


0.015 


0.015 


0.019 


0.006 


0.010 


NI 


1068.1 


0.699 


5.322 


5.297 


5.294 


5.318 


5.314 


0.013 


0.027 


0.030 


0.022 


0.026 


DI 


1068.1 


0.731 


5.343 


5.319 


5.315 


5.340 


5.336 


0.013 


0.027 


0.030 


0.022 


0.026 



Col. (1) spectrum identification: DI: Delbouille intensity, NI: Neckel intensity, NF: Neckel flux, KF: Kurucz flux. Col. (2) wavelength of the 
line. Col. (3) measured equivalent width. Col. (4) phosphorus abundance, A(P), according to the CO^BOLD 3D model. Cols. (5)-(8) provide 
A(P) from ID models, odd numbered cols, correspond to a micro-turbulence ^ of 1.0km s"' , even numbered cols, to f = 1.5 km s"' . Col. (9) 
uncertainty in A(P) due to the uncertainty in the measured EWs. Col. (10)-(13) provide 3D corrections, even numbered cols, for ^ = 1.0 km s"', 
and odd numbered cols, for 1.5 km s"', respectively. 



